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Abstract 

A systematic construction of Stackel systems in separated coordinates 
and its relation to bi-Hamiltonian formalism are considered. A general 
form of related hydrodynamic systems, integrable by the Hamilton- Jacobi 
method, is derived. One Casimir bi-Hamiltonian case is studed in details 
and in this case, a systematic construction of related hydrodynamic sys- 
tems in arbitrary coordinates is presented, using the cofactor method and 
soliton symmetry constraints. 
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1 Introduction 

There is a quite well developed theory of the passage from an integrable, infinite 
dimensional Hamiltonian system (soliton system) to its various constrained flows 
which are themselves completely integrable Hamiltonian systems. Actually, by 
using the Hamilton-Jacobi method with respect to two evolution parameters 
X and TV-gap solutions and A^-soliton solutions of a given PDE can be con- 
structed directly from solutions of related ODEs (constrained flows) 

In the present paper we are interested in, instead of soliton systems, the first 
order quasi-linear PDEs of the form 

n 

qu =^Wij{q)qj^, q^ = q,{x,t), i = l,...,n, (1.1) 
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called hydrodynamic or dispersionless systems. More precisely, we consider these 



systems among (1.1), whose general solutions can be obtained from the solu- 
tions of the related integrable finite dimensional Hamiltonian systems, like in 
the case of soliton systems. Such an approach was presented for the first time 
by Ferapontov and Fordy in , where authors gave general solutions of appro- 
priate hydrodynamic systems from general solutions of related separable finite 
dimensional systems. The idea is the following. Consider for example a com- 
pletely integrable Hamiltonian system of two degrees of freedom, given by a 
Hamiltonian function 

H=^{pl+pl) + V{q,,q2) (1.2) 
and an additional constant of motion 

F = q2PiP2 - qipj + W{qi,q2), (1.3) 

which commutes with H with respect to the canonical Poisson bracket. Note 
that both functions are quadratic in momenta, which belongs to the class to be 
considered in this paper. Let x be an evolution parameter of the flow generated 
by H and t be an evolution parameter of the flow generated by F. The commu- 
tativity of the two flows means that we can consider a 2-dimensional surface in 
phase space, parametrized by x and t. On this surface, the equations of motion 
for q are 

dH dF 

= qit = ^, 1 = 1,2. (1.4) 

OPi Opi 

Eliminating pi and p2 from ( |l.4| ), we obtain a system of hydrodynamic type 
(0) 

qit = q2q2x, q2t = q2qix ~ '2,qiq2x. (1.5) 

In this calculation, V and W play no role, so in fact the hydrodynamic system is 
generated by the geodesic parts of both Hamiltonian functions. The functions 
V, W for which H and F commute belong to the Stackel class of parabolic 
coordinates and arc the conserved density and flux for the hydrodynamic system 
respectively, as Vt = W^- Moreover, as there is an infinite hierarchy of separable 
potentials V, W, we have an infinite hierarchy of conserved densities and related 



fluxes. Hence, solutions of Hamiltonian systems (1.2) with arbitrary separable 



potential V are simultaneous solutions of hydrodynamic system (1.5) 



Remark 1 One can pass to higher order PDEs eliminating one of the qi but 
it is necessary to specify a particular form of V . The Henon-Heiles potential 

V = qf + ^(7iq| leads to the KdV equation for qi J^, whilst the quartic potential 

V = IGqf + 12qlq2 + g| generates 

1 / qixx , 2 

qit = -7^ 1- 64(ji 



48 V qi 

which is not integrable PDE ( all the integrable cases of such equations are listed 
in canonical form in 1^). 
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Here following this line we consider the problem more systematically. First, 
we observe that all examples from belong to subclass of Stackel systems 
considered by Benenti which, as was shown by Ibort, Magri and Marmo 
[ p^S] , are one-Casimir bi-Hamiltonian systems. Hence, applying one-Casimir 
and developing multi-Casimir bi-Hamiltonian formalism to quadratic in mo- 
menta separable systems, we were able to construct in very systematic way 
related hydrodynamic systems together with their general solutions. Actually, 
in Section 2 we derive explicitly, in separated coordinates, the general form of 
Stackel systems on Riemann (pseudo-Riemann) manifold and the related hy- 
drodynamic systems integrable by Hamilton-Jacobi method. Then, in Section 
3, we connect the considered Stackel systems with a bi-Hamiltonian formalism 
generalizing one-Casimir case onto multi-Casimir one. In Section 4 we present 
a one-Casimir bi-Hamiltonian systems in arbitrary coordinates (not necessary 
canonical), which we call the cofactor Stackel systems, and related cofactor 
hydrodynamic counterparts. Then, we present a recipe for the construction of 
some class of cofactor hydrodynamic systems in the Cartesian coordinate frame. 
Finally, in Section 5, we use constrained flows of soliton systems for a systematic 
derivation of hydrodynamic systems from the class considered. 

2 From Stackel Hamiltonians to complete inte- 
gral of related hydrodynamic systems 

All examples discussed in this paper belong to the class of separable systems 
associated with Stackel matrices. Actually, in 1893 Stackel gave the first charac- 
terization of the Riemann (pseudo-Riemann) manifold {Q, g) on which the equa- 
tions of geodesic motion can be solved by separation of variables. He proved that 
if in a system of orthogonal coordinates (A, n) there exists a non-singular matrix 
(p — ('y3^(Afc)), called a Stackel matrix such that the geodesic Hamiltonians 
are of the form 

n 

Er ^ Y.^^-')>1 (2.1) 

1=1 

then Er are functionally independent, pairwise commute with respect to the 
canonical Poisson bracket and the Hamilton-Jacobi equation associated to Ei 
is separable. 

Then, Eisenhart gave a coordinate-free representation for Stackel geodesic 
motion introducing special family of Killing tensors. He proved |^ that the 



geodesic Hamiltonians can be transformed into a Stackel form (2.1) if the con- 
travariant metric tensor G = g~^ has {n—1) commuting independent contravari- 
ant Killing tensors Ar of a second order such that 

admitting a common system of closed eigenforms ai 

{Al-vlG)a,^Q, da,=Q, i = l,...,n, (2.3) 
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where vl are eigenvalues of (1, 1) Killing tensor Kr = Arg {K* — gA*). In local 
coordinates g on Q we have 



Kr = Y,iKr))-^®dqj, K;^Y.^Kr)ldq,®^. (2.4) 

From now on, separated canonical coordinates will be denoted by (A, /i) and 
natural coordinates, not necessarily canonical, by {q,p). For n degrees of free- 
dom, let us consider n Stackel Hamiltonian functions in separated coordinates 
in the following form 

n 

Hr = J2vl.g"fl^,+Vr{X)^ fi^KrGfl + Vr{X), r=l,...,7i, (2.5) 

1=1 

where = {fii, n^)^ ,a.nd K(A) are appropriate potentials separable in (A,/x) 
coordinates. So, we have an n-dimensional surface parametrized by n Hamil- 
tonian "times" ti = x,t2, ■■■,tn- In this case we have ^n{n — 1) hydrodynamic 
systems (^^) written down in the form 

=<iArt,, wli^vl/vl k>l = l,2,...,n-l, (2.6) 

and the corresponding conservation laws: Vit,. — Vkti- For a given evolution 
parameter tk and a "space" coordinate t;, there are n — 2 hydrodynamic flows 
Ait,. = w^iXitnT k ^ I which commute with a given flow. It follows directly 



from the involutivity of (2.5). The requirement that H = iJi, and are in 



involution with respect to the canonical Poisson bracket leads to the equations: 
div\ = Q for any i = 1, n, (2.7) 

d, ln(g") = for any z ^ j, (2.8) 

^^Vk - vld.Vi = for any z = 1, n, (2.9) 
where di :— d/dqi. Condition ( |2.7D means that our system (2.6) is linearly 



degenerate. Cross differentiation of ( p.8[ ) gives 

(4^) = (-^) ^^j^i^^- (2-10) 

This last condition is called the " semi-Hamiltonian" property in the context 
of systems of hydrodynamic type. In that sense we consider weakly nonlinear 
semi-Hamiltonian (WNSH) hydrodynamic systems. 

The solution of the system of equations (^.7])-( pj| ) is easy to derive if one 
realizes that the separated coordinates (A, /i) for a Liouvillc integrable system 
have to fulfil the Sklyanin conditions 

•f^{X^,^ii■,Hl,...,Hn) ^Q, i = l,...,n, (2.11) 
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which guarantees the solvabiUty of an appropriate Hamilton- Jacobi equation. 
For the integrable system (2.5), under the assumption that all functions ip^ are 
linear with respect to all Hj, conditions (2.11) take the general form 



/.(A.)m? +7.(A.) = Y.'^'i{Xi)Hk, i = 1,. 



(2.12) 



k=l 



where fi,"f^,^^ are arbitrary smooth functions of its argument and the nor- 
malization <i>" = 1, i — l,...,nis assumed. To get the explicit form of 
Hk = Hk{X, fj,) one has to solve the system of linear equations (2.12). The 
results are the following 



9^' = (-1) 



,,+idetW^ 



where 



det ' 

n 

K. = ^(-l)'+'^7.(A.) 
/ $i(Ai) $2(Ai) ... 



det VF* 



(2.13) 
(2.14) 

(2.15) 

V KiXn) <^l{Xn) ••• K-HK) 1 

and W'^'' is the (n— 1) x (rt — 1) matrix obtained from W after we cancel its ith 
row and fcth column. Then the Stackel matrix ip is given by 



W = 



-(-1) 

det VF" 
detW^ ' 

<i>r\Ai) 1 



/ /i(Ai) /i(Ai) /i(Ai) /i(Ai) 



V 



^^(A^) $^(A„) 
/„(A„) /„(A„) 



■i-;;''(A„) 1 
/„(A„) /„(A„) 



(2.16) 



Notice that for r = 2, we reconstructed the result of Ferapontov for the 
functions v^. 

Remark 2 One can observe that for all known sep arabl e systems, we have fi = 
f, 7j = 7, $f = , i = 1, ...,n, so the conditions { 2.1i } are represented by the 
separation (spectral) curve 



/(A)Ai2-f7(A) = 5]a>'=(A)Hfe 



k=l 



Given a Hamiltonian system in canonical separated coordinates (A,/i) we 
can linearize the system through a canonical transformation (A,/i) — > (6, a) in 
the form bi = §§-,IJ'i = ■§^, where S{\,a) = ^^^i Si{Xi,a) is an additively 
separated generating function, satisfying the related Hamilton- Jacobi equations 



Hr{X,-^) = ar, r=l,...,n. 



(2.17) 
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For the Hamiltonian functions (2.5), fulfilling the Sklyanin conditions ( 2. 12] ), 
Si{Xi, a) are given by a system of ordinary differential equations 



/ icr \ 2 n 



fe=i 



Then, in (6, a) coordinates the flows are trivial 

(OjOt. = 0, {bj)t,. = Sjr 

and the implicit form of the trajectories Ai(ti, i„) is 



(2.18) 



(2.19) 



bj{X,a) 
where 



^l(0 



^r^(o 



= tj + const j , j — I, 



.fc=i 



1/2 



(2.20) 
(2.21) 



Multicomponent functions Xi{ti, i„) are simultaneous solutions of all dynam- 
ics defined by n Hamiltonians (2.5), as well as general solutions of all hydrody- 
namic systems ( ^.6D [ pA| . 

Of course, we would like to distinguish our separable Stackel system (2.5) 
and the related hydrodynamic systems written down in any natural coordinates 
(g,p). As we know, for one and a half century this problem has been unsolved 
in general. Quite recently, there appeared two strong formalisms which allow 
us to construct systematically a transformation to separated coordinates. One 
formalism is based on Lax representation of the system considered pC| ] and the 
other one on a bi-Hamiltonian formalism ||l^- ||l^. In what follows we will use 
some results of the second formalism. 



Bi-Hamiltonian Stackel systems in separable 
coordinates 



An important fact is that the Stackel systems (2.5) are bi-Hamiltonian on the 
phase space M — T*Q x ^ i.e. the cotangent bundle spanned by k additional 
Casimir coordinates, where 1 < k < n. Let us assume that M is equipped with a 
linear Poisson pencil IIa = Hi — Alio of rank 2n, i.e. a pair of Poisson operators 
(tensors) 11^ : T*M — > TM each of rank 2n such that their linear combination 
Hi — Alio is again a Poisson operator for any \ £ R { the operators Ho and 
Hi are then said to be compatible). Moreover, let us assume that k Casimirs 
, i = 1, k of the pencil 11^ are polynomials in A of orders ni, Uk 

i/f =i/«A"'+i7fU"'-i + ... + i7«, i^l,...,k, (3.1) 
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where ni + ... + Uk — n. By expanding equations H\{dH^^^) — 0,i — 1, ...,k 
in powers of A and comparing the coefficients of equal powers, we obtain the 
following bi-Hamiltonian chains 



no(ciijf^) 





ni(di?f ) 



: (3.2) 

= ni(di?«), 

where i = l,...,k. Notice that each chain starts with a Casimir of the first 
Poisson operator and terminates with a Casimir of the second Poisson operator. 
As follows from (3.2), the functions are in involution with respect to both 

Poisson structures. If additionally all Hj are functionally independent, then 
the chains define a Liouville integrable system. Let us introduce the following 
Casimir coordinates Ci — H'^\i — l,...,fc. From the definition, a Darboux- 
Nijenhuis (DN) separated coordinates (A,/i,c) are the canonical coordinates in 
which both Poisson structures take the following form 



n 



— 



/ 

-/ 





\ 




V 



(3.3) 



/ 



Hi = 



-A 



A 




V 



^H[^'>/^^l \ 

-dH^^/d\ 



J 



where / is n x n unit matrix, A = diag{Xi, A„) and the symbol * denotes the 
elements that make the matrix skew-symmetric. 

Now, n Sklyanin conditions ( ^.11 ), with n functions Hj^\i — l,...,k, j = 



/.(A0m^+7.(A.)=E*'(^*)^a 



i = 1, ...,n, 



(3.4) 



where fi,j^,"^l are arbitrary smooth functions of its argument, H^^ are the 



Casimir polynomials (3.1) evaluated in A = A 



(3.5) 
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and the normalization 5"*^ — 1, i = l,...,nis assumed. One can immedi- 
ately reconstruct conditions (2.12), as Hi — h[^\ H2 — ...^ Hn — Hn^J , 



n-l 



(Aj) = A,, 7,(A,) 



$KAi) = Ar-^vi;i(A,), $f(A.) = K^-'^}{Xi),..., $ 

7.(A») - Ej=i Ar'*KA^)cj and then, from (|2l|) we get H^/'> = H^'\X,tx,c). 

The simplest case of one Casimir is determined by the following Sklyanin 
conditions: 



f^{X^)^l! + 7.(A.) = cA^ + -ffiAr ' + - + H,,, z = 1, , 



(3.6) 



Here we put for simplicity ci = c, i/^. = i?fc. Notice that related hydrodynamic 
systems, written in A coordinates, are completely described by functions 
(2.13) which are determined by the right hand side of the Sklyanin conditions 
(2.12). Because in one Casimir case the r.h.s. of eq. (3.6) is fixed, there is a 
unique set of functions u^.. Actually, as was found in |12|, 



E 

i=l 



dPr /»(A.,)Atf +7,(A,) 



+ cp,,(A), r = l,...,n. 



(3.7) 



where — JJj-^j^iXi — Aj), and ^^(A) are coefficients of the characteristic poly- 
nomial of A 



det(A/ - A) = (A - Ai)(A - A2)...(A - A„) = ^ 



i=0 



i.e. are Viete polynomials. In the notation of eq.(|2.5|) it means that 
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9A.-' 2^ 



dX, A, ■ 



(3.8) 



(3.9) 



so Kr = diag{v}., ...,v"), G = diag{g^^ 
If takes the form 



/ 



if 



g"") and the related Stackel matrix 

Ai 1 \ 

/i(Ai) ■■" /i(Ai) 77(177 



A"- 



A?- 



A"" 



A"" 



(3.10) 



(A„) /„(A„) 



/„(A„) /„(A„) 



An additional term cp^(A) is related to the new Casimir coordinate and can be 
absorbed by the potential. Note that the Killing tensors Kr are given by the 
cofactor representation 



3/(A/-A)-^i^„_,A\ 



(3.11) 



where cof{A) stands for the matrix of cofactors, so that cof{A)A = {detA)I. 
The cofactor nature of Kr gives immediately the following relation 



i^.+i = ^p,(A)A'-'= 



(3.12) 



k=0 
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and vice versa, Kr given by the relation (3.12) are of cofactor form (3.11). The 
functions 7?^ ( ^/^ ) form a single bi-Hamiltonian chain (3.2), where 

/ \ /OA dHi/dfi 

Ho = ( -/ , Hi = -A -dHi/dX | . (3.13) 
000/ y*.^ 

Observe that both Poisson structures can be projected onto T*Q 
and 

^ = ^i^o'=(o I) (3.15) 

is a (1, 1) tensor on T*Q with a vanishing Nijenhuis torsion. The operator N is 
just a lift from Q to T*Q of a (1, 1) tensor A 

d 

^ = y" >^^^(^dX^ (3.16) 

on Q with a vanishing Nijenhuis torsion jl^. Moreover, as 

A*dX, = X^dX,, (3.17) 

then 

r r 

if;+idA. = ^p,(A)(A*)'^-'=dA. = ^p,(A)Ar'dA. 

fc=0 fc=0 

= = f;+idA„ (3.18) 

and multiplying both sides of eq.(3.18) by G we get 

{GK* - vi.G)dX., = 0^{A;- vl.G)dX^ = 0, i = 1, ...,n, (3.19) 

i.e. the tensorial Eiscnhart realization ( p7^ ) of the Stackel results. Moreover, if 
we define 

A := i^^AG/i, (3.20) 

then 

{A,Ei}g„ = aEi, a = Gd{Tr K) I dX. (3.21) 

A of such property is called conformal Killing tensor with the associated poten- 
tial in the form of Tr A. Finally, note that Hi and i7„ are related by 

-p„no(di?i) = ni(rfij„), (3.22) 
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which is known as the quasi-bi-Hamiltonian representation and is just a result 
of the projection of one Casimir Poisson pencil onto a symplectic leaf of XIq , 

Thus, within the class of one Casimir Stackel systems, i.e. cofactor Stdckel 
systems, infinitely many systems from that class are related toin(n — 1) hydro- 
dynamic systems (2.6) governed by n Killing matrices Kr from geodesic 
Hamiltonians. For example: v\ = l,vl — K~J2k=i and = (—1)" (Ilfc^i ^k) 

There exists a sequence of generic separable potentials Vr'''^ , k = ±1, ±2, 
which can be added to geodesic Hamiltonians, given by the following recursion 
relation Q 

V,^''+'^ ^-V}1{ + V}'Wl''\ V}'^^p,, A; =1,2,..., (3.23) 
and its inverse 

K^-''-'^ - -V}:^^ + Vy^^Vf^K K'-i) = p,_i/p„, fc = 1, 2, ... . (3.24) 

Using the notation V\ — X]j=o ^n-j^'' , the recursion formulas (|^ and ( ^^ 
can be written in a compact form 

+1) = det(A/ - A)VI''^ - AFf ^ (3.25) 

and 

Potentials F^'"'-' ( 3.23 ) and T/^"*^) ( |3.24| ) are gener ated by the corresponding 
monomials 7i(Ai) — A"^*^"^ and 7i(Ai) A^'^ from (3.6). The infinite hierarchy 

of conservation laws takes the form (^v}''^^ = (Vs^j , k = ±1, ±2, ... . 

Now, the complete integral ( 2.2Cl| ) of hydrodynamic systems (|2.6| ) with = 
-§17 is given by 

^—-d^ + ...+ ^^d^^t.+constj, j = l,...,n, (3.27) 

where the functions V3j(f) are arbitrary. Additionally, one can always construct 
a (1 + (n — l))-dimensional hydrodynamic system with the solution of the form 
( ^.27 ) involving all independent variables tj,j — 1, ...,rt simultaneously. 

Example 1 Three-field system (Ai, A2,A3). 

In this case, we have 

Pi = — Ai — A2 — A3, 

P2 - A1A2 + A1A3 + A2A3, (3.28) 

P3 = ^AiA2A3, 

hence the following three WNSH hydrodynamic systems in Riemann invariant 
form are admissible: 

Ait2 = — (A2 -l- A3)Ait^, A2t2 = "(Ai + A3)A2ti, A3t2 = — (Ai -|- A2)A3t^, (3.29) 
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^1*3 — ^2^3^1ti, ^2t3 — ^1^3^2ti, ^3t3 — AiA2A3t^, (3.30) 



'^"a ~ ~ \ ^1 ^\ ^1*2) A2t3 — - . ^ \ A2t2, A3t3 - - , \ ! A3t2 . (3.31) 
^2 + A3 Ai + A3 Ai + A2 

The complete integral for all these systems is given by 

/ ^77T'^C+ / ^T-rrd^ = +constj, j = 1,2,3. (3.32) 



In each of the above cases, only a pair of ti coordinates is involved, so the third 
one can be put equal to zero. Notice that the first two equations can be coupled 
into a (1 + 2)-dimcnsional hydrodynamic system 

1 2 2 1 
Ait3 = -2'-'^2 + A3)Aiti - -(A2 + A3)Ait2, 

1 2 2 1 
A2t3 = -2^^'>-^^^^^^*^^2^^^^^^^^^*^' 
1 2 2 1 

A3t3 = - ^'^•^l + •^2)A3ti - 2('^l + ^2)A3t2, 



for which (3.32) is an integral involving simultaneously a triple of independent 
coordinates U, i = 1,2, 3. 

The two-Casimir case is given by the following Sklyanin conditions 

MX,)^J^ + 7,(A.) = *i(ciAr + H['h-^-' + ... + <)) (3.33) 
+C2X'l'+H[^hl''-^ + ... + Hil\ 

where i = 1, ...,n and ni + n2 = n. Because are determined by the r.h.s. 
of (3.31), involving arbitrary functions we have infinitely many sets of v^. 
functions and infinitely many different WNSH hydrodynamic systems written in 
A coordinates, integrable by the Hamilton- Jacobi method. A particular example, 
with ni = 1, 712 = n — 1 and = A", can be found in jl^ . 

As mentioned before, we are interested in constructing hydrodynamic sys- 
tems and a hierarchy of conservation laws, written down in some natural coor- 
dinates, which are integrable by the Hamilton- Jacobi method. Then, we would 
like to find an appropriate transformation to Riemann invariant form (|]^) rep- 
resented by separated coordinates. In the next section, we present some results, 
based on the known theory on the separable one Casimir bi-Hamiltonian chains 
fF9|-|2l. 



4 Separable cofactor systems in arbitrary coor- 
dinates 

As most relations presented in the previous section, although derived in sepa- 
rated coordinates, are of tensorial form, i.e. coordinate-free form so they are 
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valid in arbitrary coordinate frame spanning Q. First, let us restrict our consid- 
erations to a class of one-Casimir bi-Hamiltonian systems written in arbitrary 
canonical coordinates governed by a nondegenerate point transformation be- 
tween A and q coordinates. After the transformation, one gets the following 
Poisson structures 

/ / \ / L{q) dHi/dp \ 

Ho = -/ , Hi = -L^(g) F{q,p) -dHjdq , (4.1) 
\000/ V* * 0/ 

where Fij = -^{Lp)j — -^{p^L)i and a Casimir of the pencil Ha = Hi — Alio 
of the following form 

n 

Hx{q,p,c) = ^i/„_,(g,p,c)A'^ (4.2) 

1=0 

= p'^cofilX - L)G{q)p + Vl^^\q) + cdct(A/ - L), 
where separable potentials are generated by the recursions 

+^)(g) - det(A/ - L)V^''\q) - AFf '(g), (4.3) 



vt'-'\^ - \ ( '^^[,'J^^V ^^)(,) - Vr\q)) . (4.4) 
The geodesic Hamiltonians are as follows 



n-l 

Er{q,p) = P^AMp = p^ Kr{q)G{q)p, ^ Kn-^{q)\' = cof{\I - L{q)), 

i=0 

(4.5) 

with a contravariant metric tensor G{q) and 

Kr = j2P>^il)L'^~'^ (4-6) 

fc=0 

where /0^(g) are coefficients of a characteristic polynomial of L{q). Note, that 
cofactor Stackel systems are exactly those considered by Benenti |p7| , ||l8|| , [Tstl . 
Then, the hydrodynamic equations in q representation take the form 

qt^ ^ K,K-'qt^, J > 1^1,..., n-l. (4.7) 

This class of hydrodynamic systems will be called cofactor hydrodynamic sys- 
tems. Observe, that to get cofactor hydrodynamic systems (4.7), we only need 
a (1,1) conformal Killing tensor L{q). Moreover, the separated coordinates 
Xi,i = 1, ...,n, in which the system (4.7) takes the WNSH form (2.6), are given 
by the roots of 

det(A/-L) = 0. (4.8) 
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Unfortunately, there is no systematic method of constructing L tensors in a 
general case. Nevertheless, recently some progress has been made in a special 
case of flat spaces and Cartesian q coordinates. 

Let us consider a class of one-Casimir bi-Hamiltonian systems on a flat space 
Q = i?", introduced recently in Here we briefly review the results 

which are important for our construction. Let q — {qi, ...,qn)'^ be a set of 
Cartesian coordinates and A, annxn matrix, whose elements fulfil the following 
equations: 

^^AJk + dgAu + dkA,.j = 0, i, j, fc = 1, n. (4.9) 



The equations (4.9) imply that the matrix A is a Killing matrix. An important 
class of solutions of these equations have the form |^ 

A = coJ{G) (4.10) 

with 

G = aqq^ + Pq^ + qf + (4.11) 

where a is a real constant, /3 = (/3;^, /3„)-^ is a column vector of constants and 
where 7 is a symmetric n x n matrix. One can show that for n = 2 it is the 
general solution of (4.9) 
Now, let 

G ^aqq'^ +]3q'^ + q]]'^ + J (4.12) 



be another matrix of the form (4.11) (we assume that at least one of the matrices 
G and G is nonconstant), then, all matrices Ai,i = 1, n, defined as coefficients 
in the polynomial expansion of cof{G + AG) with respect to the real parameter 
A : 

n-l 

cof{G + XG) = Y,^n-^y (4.13) 

1=0 

with Ai = cof{G), An = cof{G), are Killing matrices. 

Let us assume that over some region of Q detG ^ 0. Then L = —GG~^, 
considered as (1, 1) tensor field on Q = i?", has a vanishing Nijenhuis torsion, 
and so it can always be diagonalized. Moreover, on T*Q we define n functions 

E,^p'^A,p, i^l,...,n (4.14) 

and a tensor 

0=(g "i?)' (4.15) 
where (pi, ...,pn)'^ are momenta coordinates and an n x n matrix F is defined 

by 

F = Np^-pN^, N^aq + f3. (4.16) 

Theorem 3 Assuming that over some region of Q the operator L has n func- 
tionally independent nonzero eigenvalues: 
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(i) Q is a Poisson tensor of rank 2n. 

(ii) Ei are functionally independent and in involution with respect to Q. 

Proof. Operator Q is skew-symmetric and the Jacobi identity can be proved 
by inspection. Moreover det G ^ guarantee its maximal rank. On the other 
hand, the cofactor form of geodesic Hamiltonians together with functional in- 
dependence of eigenvalues of L operator means that we have the Eisenhart 
representation and hence the Stackel ones in separated coordinates. ■ 

Obviously, we have Liouville integrable systems for geodetic motions, written 
in noncanonical coordinates {q,p)- Choosing G = —I we get a special case 
of canonical representations. Before we separate the system, let us construct 
related hydrodynamic systems. Equations of motion for q coordinates are as 
follows 

qu =-2GA,p, i = l,...,n, (4.17) 

hence, elimination of p coordinates, together with the relation Ai — cof{G), 
leads to the following cofactor hydrodynamic systems 

qt^ = A^'AjAr^A^qt, = GA,Ar^G-'qu, j>i^ 1, - 1. (4.18) 

In particular case ti ~ ti — x we have 

qt^ = A-^^Ajq^ = -^^GAj<lx. ] = 2, ...,n. (4.19) 



To find separated coordinates for the geodesic Hamiltonians (4.14), we have 
to put them into a bi-Hamiltonian form. It can be done on the extended phase 
space M = T*Q x R with local coordinates {q,p,c). Let us introduce functions 
Dj{q) as coefficients in the polynomial expansion of det(G + AG) 

n 

^D„_i(g)A* = det(G + AG) (4.20) 

i=0 

so that Dq — det G and D„ — det G. 
Theorem 4 Functions 

Hr^Er + c^, r = 0,...,n, Eq ^ c, (4.21) 

constitute a bi-Hamiltonian chain with respect to a pair of compatible Poisson 
structures 



n 



— 



Hi = 



(122) 



14 



which starts with the Casimir Hq = c of the first Poisson structure Hq and 
terminates with the Casimir Hn — En + Cjy^ of the second Poisson structure 
Hi. The last column o/IIi is a first vector field from the hierarchy: Hfj{dHi). 

Having a bi-Hamiltonian chain, one can systematically construct separated 
coordinates j22j. The result is given as follows: 

Theorem 5 For a geodesic Hamiltonian system ( (.14 A separated 
coordinates Xi (q) are the roots of the equation 



det(A/ - L) = ^ det(AG + G) 
and related momenta fj,^{q,p) are given by the equations 



(4.23) 



1 n'^cofiG + \{q)G)p 



'^n'^cof{G + \i{q)G)n 



, i = l, ...,n 



(4.24) 



where J7 = (N + LN) and N — aq + f3. The Poisson operators (4-22) attain the 
form 



— 



Hi 




A dHi/d^i 
-A -dHi/dX 
* * 



(4.25) 



while the Hamiltonians ( 4-21 ) have the form 

dPr fk{^k) 2 , (x\ 1 

-ffr- = -2^— A, l^k+CPr{X), r=l,...,71. 



k=l 



dXk Afe 



(4.30) 



The above theorem ensures us that under the transformation given by the 
roots of (4.23), cofactor hydrodynamic systems (4.18) turn into the WNSH form 
(2^), where {KrY^ = Sijvl = — To complete the picture from the point of 
view of Stackel systems, let us end this consideration with the following result. 

Proposition 6 The generic separable potentials of the cofactor Stackel system 
(4-21), (4-22) are given by the recursion formulas |^/ 



(fc+i) 



^^^^^±^vl''\q) - XVt\,) 
det{XI - L)vl''\q) ^ Xvl''\q), 



(4.31) 
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(-fe-1) 



1 / det{XG + G) 



A 



detG 



(4.32) 



^ 1 / dct(A/-L) 

A V dctL " ^ 

where V\ —T^^q Vn-i'^^ md the potential Vr belongs to the appropriate Hamil- 
tonian Er ( 4-21 )■ 

In the canonical coordinates G = —I, G = L and we reconstruct the poten- 
tials (gj) and (lU). 

Example 2 A two field cofactor hydrodynamic system is defined by 



G 



We have 



1 91 
2^2 -qi 

-qi 1 



A. 



ql + i qiq2 

9192 ?! 



9i 



-gi92 



and hence, form (4.19) 

'72*2 = 



92(92 - g?) 



gi'?2 9? + 1 

91(92 - 9i - 1) 



2q2-qt 2^2-9? 
92(91+2) , 92(91+2) 

2-1Ui + -TT 72-9lti 



91*1, 



292 - qf 



292 - qf 



The transformation ( [1.23D to the WNSH form 

Alta = A2Altj, A2t2 = AiA2ti, 

is given by 



91 



^v/-^i^2(Ai + l)(A2 + 2) 

^ : , 92 



A1A2 



Ai + A2 + A1A2 Ai + A2 + A1A2 

Example 3 Three field cofactor hydrodynamic systems are defined by 



G = 




G = 



We have 



293 92 9i 
92 -1 
qi -1 



91 
-92 




92 -293 



-92 
293 + 9i92 
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and hence, the following hydrodynamic systems 

9lt2 = -91 9lti + '?2'72ti + 2(739341, 

9it2 = -9iti - 2gig2ii +9293*1, 
91*2 -9292*1 - 91 93*1, 



91*3 = 9291*1 - (91 92 + 293)92*1 + 9293*1 , 
91*3 = 91 91*1 + 9i92*i - (9192 + 2g3)g3t,, 
91*3 = 91*1 - 9192*1 + 9293*1, 



91*3 



-91*2 



1 2g| + 3qlq2 + 2gi g3 

2 93 - 91 92 - 9? 
^ 1 4919293 + 2gi + 

2 93 - 91 92 - qf 



1 29293 + 29^93 + 9^92 



93 - 91 92 - 9i 



-92*2 



9i 



9?9i 



-93*2, 



The transformation ([4.231) to the WNSH form ( p.29D -( p.3lD is given by 

1 



9i 

92 

93 



1 



1 



7^2 



pI 



7P1P2 



16 



pI 



where Pj are define by (3.26). 



5 Hydrodynamic equations related to constrained 
flows of soliton systems 

As it was mentioned in the introduction, symmetry constraints of soliton systems 
give us a systematic method of constructing Liouville integrable finite dimen- 
sional Hamiltonian systems Moreover, most of the examples constructed 
so far, belongs to separable systems, separated either through a bi-Hamiltonian 
formalism Q , |^ or through a spectral curve method Q , ||] . If additionally in- 
volutive constants of motion are quadratic forms of momenta, then we can again 
systematically construct related hydrodynamic systems which have a complete 



integral in the form (3.27) 



We illustrate the approach by two examples. The first one is related to 
the constrained Schrodinger spectral problem of the KdV hierarchy, where the 
natural coordinates are its n eigenfunctions. Let q^ be an eigenfunction of the 
KdV Lax operator + u with an eigenvalue ak 

qkxx +uqk ^ akqk, fc=:l,...,n. (5.1) 
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Under the symmetry constraint 



1=1 



E 



(5.2) 



of the KdV equation, where c plays a role of the additional Casimir coordinate, 
we obtain the Garnier system, well known in the classical mechanics: 



qk 



(5.3) 



The bi-Hamiltonian representation of (5^), in canonical coordinates, was found 
in M, [E5[, where the second Poisson structure is 







A 



1 T 

^ /-p A ^ T 1 T 



P 

[A^cl- {q, q)]q 




(5.4) 



Pk = Qkx, fc = 1, ...,n, A = diag{ai, ...,q;„) and (g,^) = q'^q = X]"=i Observe 
that this is again a cofactor Stackel type system where G — —I and G = 
L = A — ^qq^ ■ Geodesic motion is separable in generalized elliptic coordinates 
Ai, A„, defined by the relation 

n;u(^-A,) 



fc=i " 



(5.5) 



which are just the DN coordinates defined by the roots of det(AG' + G) — 
[ p!2| . The Garnier potential is the first nontrivial one from generic potentials 
generated by the recursion (4.3), all are separable in elliptic coordinates (5.5). 
The geodesic Hamiltonians are of the form 

k 

Hk+1 - \{P.A^P) + \Y,[{q,A=-\){p,A^-^p) - {q,A^-^p){q,A^-=p)l 
J=l 

where A; = 0, ...,n — 1 and they allow us to construct related hydrodynamic 
systems. For n = 3 the three cofactor hydrodynamic systems are 

Mq-x, qt3 = ^39a:, 9*3 = ^3^2 ^%2: (5-6) 



where 

A2 = 



qt2 



\qk + 593 ~0i2-0i-i 

-\q\q2 
-^q'^q?, 

/ 19 1? 

' -2"3g2 - 2^293 

+Q;2a3 



A, 



^(^■iqxqi 

\oL2q\q-i 



~\q\q2 

hi + hi - ai - "3 

-59293 



2"39l92 

19 12 
-2*^391 - 2"l93 

+aia3 

5^19293 



-|9i93 
-59293 
59i + 592 - "1 



2a29i93 

5ai9293 

-5a29? - 5019! 
+Q!ia2 



a2 
\ 
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The system ( ^.6[ ) takes a WNSH form ( 3.29| )-(3.31) in the generahzes elhptic 
coordinates after the following coordinate transformation 

Pi(A) = -^{ql + ql + ql) + ai + a2 + a3, 

P2W = ^[(Q!2 + "3)9? + (ai + Q!3)'72 + ("1 + "2)93] - (q:iQ;2 + aias + Q;2a3), 
P3W = -^{c^2a3qf + aia^ql + aia2ql) + aia2a3. 

Of course just the KdV case offers us infinitely many Liouville integrable bi- 
Hamiltonian systems, generated from the symmetry constraint of the KdV equa- 
tion 

qkxx + uqk = auqk, k^l,...,n, 

-T— = + y^qi + c, (5.7) 

ou ^-^ ou ^-^ 

k=l k=l 

where hm is the m-th conserved functional of the KdV hierarchy. In each case, 
we get a cofactor type hydrodynamic system. 

Just to demonstrate a vast universality of this approach, our second exam- 
ple is related to another soliton hierarchy, represented by the Jaulent Miodek 
spectral problem ||2^ (a special case of Antonowicz and Fordy spectral problem 

[0) 

Pi i V af-uia^-UQ J \ 



i = 1, n, (5.8) 



and its m-th symmetry constraint 



Shm _ 1 / iq,Aq) + c 
6u ~ 2\ {q,q) 

Here we consider the case of m = 4 [| 

7 r 5 5 3 1 

^ 128"^ ^ 16"^"° ~ 32"^"^"^ ^ 8^°"^ ~ gWla:MO:r- 

By introducing the Ostrogradsky coordinates 

<7n+l = Ul, qn+2 = Uo, 

SHi 5 1 SHi 1 

ouix lb 8 ouqx 8 



(5.9) 



Eqs. (5.8), ( |5.9D for m — A are transformed into a canonical finite dimensional 
Hamiltonian system with the Hamiltonian function 

Hi = ^{p,p)-^{q,Aq) + ^qn+i{q,Aq) + ^qn+2{q,q)+cqn+i 

-&Pn+lPn+2 + ^Qqn+lPn+2 ^ Y^'ln+l<ln+2 ~ g'?Ji+l9n+2 ~ J^^""*""^' 



19 



The bi-Hamiltonian (quasi-bi-Hamiltonian) representation was found in [|, H 
with the second Poisson structure Ho (8o) in the form ( [4.1| ) with L{q) matrix 
of the size (n + 2) x (n + 2) 



L 



A 

Olxn 



q-n+i 1 



which gives us a set of (n + 2) component cofactor hydrodynamic systems. Let 
us look at the three-component case of n = 1, then 



92 

1 n 15 „2 




1 

192 



and cof{XI - L) = /A + A2A + ^3, where 

1 

/ 2^2 

A2- 




111 



2<l2 - a 

A 

8 



292 -^92 - 593 -92 - a 



A3 = 



hi + 193 
292 
-29i 



n T"'?2 - 59192 + ^aq3 aq2 
and the related cofactor hydrodynamic systems are 



19192 
~|a92 



-491 
—a 



9*2 



^29a:, 9*3 = ^39x, 9*3 = -43^2 ^9*3 



(5.10) 



The system (|5.10D takes a WNSH form (|3.29D - (|3.3lD after the following coordi- 
nate transformation 

91 = [^a^-5aVi(A) + ^ap?(A)-ap2(A)-p3(A)]/[2a-2pi(A)], 

92 = 2a-2pi(A), 

93 - -a2 + 2api(A)-3p}(A)-2p2(A). 



6 Concluding remarks 

In this paper, developing the ideas of Ferapontov and Fordy Q and Ibort, Magri 
and Marmo jl^, we have established the relation between the bi-Hamiltonian 
family of Stackel systems and the class of hydrodynamic systems whose complete 
integral is constructed from a complete solution of the related Stackel family. 
Moreover, we have found the most general admissible Riemann invariant form 
of such hydrodynamic systems in separated coordinates. In a particular case 
of one Casimir Stackel family (cofactor Stackel systems), we also presented 
some systematic methods of construction of related hydrodynamic systems in 
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arbitrary coordinates and the explicit form of the transformation to the Riemann 
invariant form. The second of the presented methods reveals a new interesting 
link between soliton systems and the class of hydrodynamic systems considered. 
Actually, the link is as follows: 

soliton system Stackel system — > hydrodynamic system 

and the complete solution to the Stackel system is a simultaneous particular 
solution of both soliton and hydrodynamic system. 
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